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ABSTRACT 

^e have observed pulse phase-depernWnt spectral changes in the 
high energy (>2G keV) continuum of Hercules X-1. Cyclotron absorption 

:Y 

ji 

of underlying continua ir>an reprodluce the observed angular dependence 
in the high energy cutoff. We discuss implications of this model ,_which 
include the possibility of determining the angular separation l|^tween 
the line of sight a-:v} the neutron star magnetic field if the a_^ 5 orbing 
electron spectrum is known. 
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I. INTRODUCTION 


Hercules X-1 is an X-ray pulsar in a binary stellar system {Tananbaum 
et aj^. 1972), It exhibits a variety of temporal phenomena: 1.2- second 

pulsations, 1,7-day binary period, 35-day high-low cycle (cf. Giacconi 
1973; Davison and Fabian 1976) — and spectral features (cf. Clark 
1972; Ulmer 1973; Holt ^ al,. 1974; Becker e^ 1977; 

Pravdo ^ 1977a, Paper 1; Pravdo et aj.. 1977b, Paper 2). This communi- 

cation reports on observations of Her X-1 performed by the Goddard Space 
Flight Center Cosmic X-ray SpectfOs copy Expert (GXS) on QSQ-& during , 



August 31 - September 6, 1977. This was the second observation of the 
source with the CXS. A detailed study of the 2-20 keV spectrum of Her X-1 
as a function of 1.2 second.pulse phase was among the results (Paper 1) 
of the first observation in August 1975. Paper 1 reported the existence 
of a small-region 0.16) of pulse phase, asymmetrically located within 
the pulse peaks, in^which the spectrum progressly hardens and then softens 
to its initial off-lpulse value. We have extended this investigation by the 


use of a high energy CXS detector which is sensitive to X-roys^between 2 
and 60 keV. We report here pulse-phased spectroscopy > 20 keV, and demon- 
strate the consistency of all the 2-60 keV data with a cyclotron-absorbed 
X-ray source. 

II. EXPERIMENT AND ANALYSIS 

\V \\ 

The detector is a multi-layer multi -anode xenon-filled proportional 
chamber with a mechanically collimated 5° FWHM circular field of view 

i ^ '> 

(cf. Serlemitsos e _ t al. 1976). It is pointed 5° off the negative spin 
axis so that a 5° cone is swept out during each satellite spin ("“ 10 sec). 
The scan enables us to obtain a nearly simultaneous observation of the 
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diffuse X-ray background in theff^ource free regions ^jacent to Her X-1. 
The net Her X-1 spectrum is determined after subtracting off-source data 
from on-source data. Each pulse-height-analyzed (PHA) count is timed to 
within 20 msec. Thus an individual PHA count can be placed in itr^APpro- 


priate pulse phase bin after making time-of-flight corrections for the 


0 


Her X-1 binary orbit, the satellite orbit, and the Earth's motion. The 
pulse period of 1.2377967 + 0.0000004s {Pravdo^^. 1977c) has been 
divided into 62 temporal bins. 


The 


of spectral analysis has been presented elsewhere 


(Serlemitsos ^ al_. 1975). \n brief, an analytical model of the incident 


specvtrum is multiplied by a previously determined detector response matrix. 




This resultant is then compared to the net PHA data using a x test. 

Once an acceptable fit is found, an inferred incident spectrum is obtained 
by dividing the PHA data by spectrum-dependent channel -by-“Channel effic- 


iencies. In' pji^)ctice we have used a computer automated modification of 




the CURFIT subroutine described in Bevington (1969). This program yields 
estimates of errors on fitting parameters, which are close to one sigma 
errors whenever the "curvature" search option is used (Marquardt 1963). 

Unless otherwise stated these one sigma errors are presented. In other 
cases, the criteria discussed by Lampton et (1976) and Avni (1976) 

p 

for increases in a minimum value of x are employed to determine parameter 
bounds with higher than one sigma confidence. 

We have employed tv/o basic models of the incident spectra for analysis. 
The first is a simple empirical moder containing four parameter^. It is 
of the form: 
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^ (photons cnr^sec~^keV E for E 5 , P 3 

E-P 

= exp-{-^) for E > P 3 . 

~r 

This model serves in presenting the general t»'«jnds of the spectra, 
but has no physical interpretation. The second model also contains 
an underlying power law spectrum with a high energy cutoff due 
to cyclotron absorption in a high (> 10^^ gauss) magnetic field (Daugherty 
and Ventura 1978). The physical parameters included in this model are 
the surfacC magnetic field strength, the angle between the line of sight 
and the magnetic field lines, the gravitational red shift, the polarization 
distribution of the X-rays, and a characteristic electron momentum, momentum 
spectrum, and column density of the absorbing electrons. Many of these 
parameters were held fixed at reasonable values (Section IV) so that 
in effect a three parameter model was fit (Section III). 

III. RESULTS 

In general there is excellent agreement with the results described 
in Paper 1. Figure 1 illustrates the Her X-1 pulse profile and spectral 
parameters in the simple empirical model as a function of pulse phase. 

The region of spectral hardening is well defined by the plot of spectral 
number index (Figure Ic). Phase zero of the pulse is defined to correspond 
to the bin in which the hardening reaches a maximum (i.e. number index 
reaches a minimum). Note that the high energy peak (> 25 keV) is asyitmetric 
with respect to phase zero, and of greater phase duration than the hardened 
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region. This strongly supports the conclusion in Paper 1 that the intensity 
changes in the pulse are at most loosely connected to the spectral changes. 

The power law index is determined and constrained largely by the 
low energy {< 15 keV) data. From the previous Her X-1 observation (Paper 
1) we determined that the < 4.5 keV continuum in off-pulse spectra indicates 
a small amount of absorption by cold material. Also, when tli)e flattened 
continue of the on-pulse spectra are extrapolated to energies below 4.5 
keV, a low energy excess results. To avoid these complications 


data below 4.5 keV are not included i n this ana lysis. 


The final two plots in Figure 1 illustrate the cutoff parameters. 

"Off “pulse" bins are defined to be all bins with plus or minus phase greater 
than 0.08 (i.e. outside the hardened region). The off-pulse spectra have a 
lower cut off energy and higher e-folding energy than the on-pulse specta. 
This fact is quantitatively shown in Figure 2. The first two plots are a 
blow-up of the on-pulse high energy parameters. Also shown are the 99% joint 
confidence limits of these parameters determined from the combined 
off-pulse spectrum. In all eases, the model fits were acceptable with 

O 

reduced chi squared (x^) 1 for 47 degrees of freedom. The net 

counts in these spectra range from 4 - 8000 for individual bins, to 
175,000 for the combined off-pulse spectrum. The spectrum averaged over 
pulse phase has ±0.3 keV and ± O.S keV 

in good agreement wibf;) the same parameters found by Becker et al . (1977). 

The final Figure shows contrasting spectra, one from phase zero, 
and the other from an edge of the on-pulse region. We find no supporting 
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evidence for line emission near 60 keV {Trumper ^ 1978; Coe ^ 

al . 1977), although our one sigma upper limit is consistent with the line 
photon intensity observed by Trumper et al. 

IV. DISCUSSION 

The lo2 second pulse period of Hercules X-1 identifies it as a 
neutron star. During stellar rotation, periodic views of an X-ray 'hot 
spot' create observed pulsations. The X-ray energy source is the 
gravitational potential energy of accreting matter. Strong stellar mag- 
netic fields cohstfain this matter to fall along field Tines to the 
magnetic poles on the surface (Lamb ^ 1973), We suggest ITiat the 

X-rays which comprise the spectrally hardened region of the pulse originate 
at the X-ray hot spot and are largely unaffected by scattering or absorption 

processes away from the hot spot (see also Paper I), in contrast to the 

/; 

X-rays observed at other pulse phases. Thus analysis of the hardened X-rays 
in this "spectrally defined. beam" represents a deep probe in toward the 
neutron star surface. 

/\ 

It is clear that significant changes occur in the < 20 keV X-ray con- 

tinuum of the above-mentioned beam. We leave discussion of thia effect 

for a future work (Bussard and Pravdo 1978). In this discussion we will 

focus on the pulse-phase dependent changes in the > 20 keV continuum. 

Boldt et a^. (1976) s|iowed that energy-dependent Compton scattering 
1 3 ’! 

in a field of ^ 10 gausscouTd reproduce the high energy cutoff. However 
the recent, lower estimates for the field stre ngth If (Trumper et al. 1978) 

!f ; 

taken together with the expected electron momenta d)i stri butiton (see below), 
indicate that cyclotron absorption dominates over Compton scattering above 
20 keV. In addition we consider the present mod'eT superior in that it 
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\V 

allows for the newly observed angular dependence in the cutoff. Cyclotron 
absorption is| the inverse and co-existant process to cyclotron emission 
(Basko and Sunyaev 1975; Daugherty and Ventura 1977; Meszaros 1978) which 
may have been observed in the Her X-1 spectrum (Trumper ^ al_. 1978). 

The X-ray emitting region is located near the surface, so that the 
magnetic field does not change much from its surface value within the 
hot spot* The electron distribution here consists of infalling electrons 
and those which have been stopped by interactions in an accretion mound 
or atmosphere (Davidson 1973). High momentum electrons will observe continuum 
photons to be D^ippler-shlfted in the electrons' rest frames. If an electron 
momentum, P satisfies,, the equation, 

c 



then the electron resonantly absorbs the photon of energy E as a cyclotron 

2 3 

line photon; where mis the electron rest mass and = 4.414 x 

13 ^ 

10 gauss. The resulting X-ray spectrum has the form 


where Pg is the number index, is the column density of electrcfts at 
the appropriate momentum ^(E,9), is the cyclotron absorption 
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1 , 0 

cross section (Oau^l[ierty and Ventura 1978) for polarization state e, 

and f . the fraction in that state. We assume that the X-rays are 

6 

unpolarized (Novick 1977). Since the electrons are contrained to 

follow field lines, the angle 0 is between the photon propagation direction 
(line of sight) and the polar field lines. 

c/ 1 2 

The magnetic field B in equation 1 must be between 2-7 x 10 gauss 

in this model. These limits are based on the observation that the high 

energy cutoff begins near 20, keV. If the magnetic field were smaller 
1 ? 

than 2 x 10 Gauss then the cutoff would begin at much lower energie.s, 



A field larger than 7 x 10 Gauss would require relativistic electrons 

•.j. 

for absorption of photons near 20 keV, We will not consider this possibility 
because other processes (e.g. pair production) become important in this 
regime. A value of 6 x 10^^ gauss is adopted for B.' This value is consistent 
with general theoretical estimates for a neutron star magnetic field (cf. 

Lamb 1974) ancj^with direct (Trumper et al. 1978) and indirect (Paper 2) 
observational estimates. A gravitational redshift of 0.1 is assumed (cf. 
iWddleditch and Nelson 1976). 

The remaining critical parameters describe the absorbing electron 

distribution. We have chosen to let the angle e be a free parameter for 

several electron distributions. This is illustrated in figure 2c and / 

2d for an exponential electron moment^um spectrum and a one-dimensional 

Maxwellian spectrum respectively. The spectra, with one exception, are 

2 

acceptably fit by the model with Xr = 0.9-1. 4 for 51 degrees of freedom. 

For the phase zero Spectrum and the Maxwellian model, Xr = 1*7, because 
the model is more sharpely cutoff than the data. In all the models tested 



the best fit values for a perform a 10° 20° excursion, with the angle 

of closest approach itientified with phase zero. The actual values of the 

angles ar^ strongly' dependent on th^-: column density of high momentum electrons. 

The low column density model (Figure 2c) could be appropriate if the absorbing 

electrons represent a small fraction of the total hot spot electrons— e.g. 

the infalling electrons. The high i lunm density model (2d) could occur 

if all the* electrons contribute to absorption. Basko and Sur\yaev (1975), 

in a model for pulse formation, estimate that energy is deposited by 

-2 

inf ailing matter within 50 gm cm .. This yields an electron column density 
similar to those employed with the 

These column densities imply that the ranges of infalling particles 
are small enough 10 cm) so that the assumption of constant field 
strength is justified. In addition, estimates of the effective 
temperature of the hot spot 30 keV) based on the X-ray continuum 
spectrum (Ulmer ^ 1972; Holt et 1974; Becker £t 1977; this 

paper) yield a similar upper limit for the plasma scale height. The ' 

characteristic electron momenta in our models were cho'k^n to be in this 
temperature range. The plots in Figure He and 2d are drawn for charact- 
eristic electron momentum = lao keV/c, and the phase zero results for 
160 keV/c and 200 keV/c are also included. : .v 

We note that the electrons which contribute to the absorption process 
have momenta as large as' A. 4 P^. These electrons on the exponential tail 
of the distributions have sufficient velocities to abss#/ ^,,20 keV photons 
in a field with a 60 keV characteristic line energ^r^ Electrons with this 
high momentum can be created in knock-on collisions with infalling ions 
(Bussard and Ramaty 1978, G amma Ra y Spectroscopy in Astrophysics , proceed- 
ings of NASA/6SFC Conference, in preparation) . 


;r » 
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Off-pulse spectra have not been specifically addressed here. There 
are no obvious changes in spectral parameters over this region of pulse 
phase (Figure 1). We find that individual off-pulse spectra are acceptably 
fit by the best-fit model of the combined off-puls^ spectrum with only 
the normalization left as a free parameter. This flux could be largely 
isotropic and originate from a Compton scattering shell of material at 
the Alfven surface (McCray and Lamb 1976; Basko and Sunyaev 1976; Ross 
al. 1978). 

This analysis can not unambiguously determine the separation between 
^e line of sight and the m^^gnetic field because the absorbing eTectron 
spectrum is unknown. A self-consistent model for the hot spot X-ray emis- 
sion ji^s needed. The gross variations in the continuum spectrui^ as a 
function of pulse phase (see Figure Ic) should be valuable in working 
backwards to the emitting electron spectrum. rAn- improvement to the present 
model would consider the effects of re-emission which can be important in 

li 

addition 'jo absorption. This work is currently in progress (Bussard 
and Pravd) 1978). 

-finally, we consider the implications of this analysis on the Her X-1 
beaming mechanism. The beam is spectrally defined and the spectra in the 
beam indicate via parameterization that the angle of closest approach to 
the magnetic pole is at phase zero. A "pencil" beam is therefore favored. 
We note that the lower energy (< 20 keV) data alone could not have antici- 
pated this model, and that all data 2-60 keV of Which we are aware may 
now be consistently reconciled. 

We thank J.L. Robinson-Saba and Dr. J.K. Daugherty for helpful 
discussions. We also thank Or. M. Lamp ton for aid in parameter space. 
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FIGURE 1 


FIGURE 2 


FIGURE 3 



FIGURE CAPTIONS 
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a. The energy-integrated pulse light curve of Hercules X-1. 

b. The high energy pulse light curve. 

c. ,d.,e. Spectral parameters in the empirical model for 

infjividual pulse bins. 

(1 

a.,b. Oh-pulse spectral parameters in the empirical model. 
The diamonds are the 99% joint confidence limits for the 
parameters of the combined off-pulse spectrum. 

c. Angle between the line of sight and the neutron star 
magnetic field versus pulse phase in the cyclotron 

absorption model. With an exponential electro, 

• ■;/' 

momentum spectrum. 

d. With a Maxwellian electron momentum spectrum. 


Incident spectra of Hercules X-1 from minimum and maximum 
(absolute value) on-pulse phase. 
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